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Abstract 

This article examines the properties of positive solutions to fully 
nonlinear systems of integral equations involving Hardy and Wolff po¬ 
tentials. The first part of the paper establishes an optimal existence 
result and a Liouville type theorem for the integral systems. Then, 
the second part examines the decay rates of positive bound states at 
infinity. In particular, a complete characterization of the asymptotic 
properties of bounded and decaying solutions is given by showing that 
such solutions vanish at infinity with two principle rates: the slow 
decay rates and the fast decay rates. In fact, the two rates can be 
fully distinguished by an integrability criterion. As an application, 
the results are shown to carry over to certain systems of quasilinear 
equations. 
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1 Introduction 

This article studies the following class of fully nonlinear systems of inte¬ 
gral equations with variable coefficients involving Hardy terms and Wolff 
potentials, 

r n(x) = ci(x)Hy^.y(|?/piu'^)(x), X € M"", 

\ v(x) = C2{x)Wp^^{\y\^'^vP){x)^ X € M”. 
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The Wolff potential of a non-negative Borel measure /r is defined by 




/^(Bt(x))\ dt 


where x € M"’, n > 3, 7 > 1, /3 > 0, /Sy < n and Bt{x) C M"’ is the open ball 
of radius t centered at x. Thus, if dfi = f dx where / € and / > 0, 

then 


WisMix) 


Jo V r-/37 ) t ■ 


Convention: Unless further specified, when considering system (jl.lj) 
we always assume the following: 


f P,q>0, 'j e ( 1 , 2 ], CTj belongs to (-/ 37 ,oo), 

I and the coefficients Cj(x) are double bounded functions. 


Here we say a function c{x) is a double bounded function if there exists a 
positive constant C such that 1/C < c{x) < C for x € M”". We say (n, v) is a 
solution of system (ll.lh if u, u € are non-negative and satisfy the 

integral equations for a.e. x € In addition, for a given positive solution, 
we say it is a decaying solution if there exist positive rates 0i and O 2 
such that u{x) ~ |x|“^i and v(x) ~ \x\~^^. Here, the notation /(x) ~ 
means there exists a positive constant c such that 1/c < |x|®/(x) < c as 
|x| —>■ 00 . 

The goals of this paper are to establish some new results concerning the 
optimal existence and non-existence of positive solutions and to continue 
the study from m on the decay properties of positive solutions for system 
Our study on the integral equations involving the Wolff potential has 
roots in the qualitative analysis of elliptic partial differential equations aris¬ 
ing from nonlinear analysis, calculus of variations, conformal geometry and 
mathematical physics. The prototypical example is the semilinear equation 
with weight, 

— Art = jxl'^u^, u > 0 in M”, (1.3) 

where p > 1 and a G M, and we illustrate soon below how this equation is 
equivalent to a very simple case of system m- Indeed, equation (II. 3p arises 
as an important stationary model for stellar cluster formation in astrophysics 
m- Liouville type theorems for (|1.3p . when combined with blow-up and 
rescaling arguments, lead to a priori estimates to a family of elliptic bound¬ 
ary value problems m, and the classification of its solutions when a = 0 
and p = {n + 2)/(n — 2) also plays an important role in the Yamabe and 
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prescribing scalar curvature problems. It turns out that the decay properties 
of solutions to equation (jl.3l) are closely related to these other properties, 
and we have a fairly good picture of the asymptotic behavior of the bound 
states. We outline such results here for completeness sake and because the 
main theorems of this paper can be viewed as generalized versions for the 
Wolff type integral systems. It is known that equation (II.3p has no solution 
if cr < — 2 or when a > —2 and I < p < (see [HI ETj). Thus, if a 
solution exists, then p > and a > —2 necessarily hold. This implies 
that < n — 2 and this motivates the following terminology for the fast 
and slow decaying solutions. It is known that bound state solutions which 
vanish at infinity must do so with two principle rates of decay: the fast rate 
u{x) ~ or the slow rate u{x) ~ |x| WT [23] (see also [Tl [TTl [T8] 

and the references therein). 

Let us make the connection between elliptic partial differential equations 
and the Wolff type integral equations more apparent. To do so, we first note 
that the Wolff potential VI(g ,|,(-) reduces to the well-known Riesz potential 
/„(•) multiplied by a positive constant when /3 = 2 and 7 = a/2. Namely, 


W^,2{f)ix) = 

2 [n — a) 


fiy) 


■dy =: 




\x — (n — a) 

System dni) under constant coefficients therefore includes a weighted vari¬ 
ant of the Hardy-Littlewood-Sobolev (HLS) type integral system 





\x-y\n-a 

\x_y\n-a 


(1.4) 


In the unweighted case, i.e., cri,cr 2 = 0, these are the Euler-Lagrange equa¬ 
tions for a functional associated with the best constant in the HLS inequality 
[24] . If p = q, ai = a 2 and u = v, the integral system reduces to an integral 
equation, which is also associated with the Euler-Lagrange equation for the 
sharp Hardy-Sobolev inequality [2l[25[[38]. The HLS type integral systems 
are naturally associated with systems of differential equations. Eor exam¬ 
ple, when a = 2A: is an even integer, ai,a 2 < 0 and p,q > 1, system dm) 
is equivalent to the poly-harmonic system of the Henon-Lane-Emden type 

[31^: 

J (—A)^u = . 

\ {—A)^v = 

which reduces to equation (11.31) if k = 1, a = ai = a 2 , p = q and u = v. 
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Remark. If k = 1, system dlSD is often ealled the Henon system when 
cri,cr 2 > 0, the Lane-Emden system when ai,a 2 = 0, or the Hardy system 
when CTi, CJ 2 < 0, but we will just refer to it as the Henon-Lane-Emden system 
in any ease. 

Similarly, if /I = 1, system dni) is closely related to the system of 7 - 
Laplace equations of the Henon-Lane-Emden type 

( —div (|Vu|'>'“^Vtt) = 

\ —div (iVup^^Vu) = C 2 (x)|x|°' 2 tt^, 

and we shall describe their close relationship in more detail shortly below. 
Other relevant examples include more general quasilinear systems, including 
those involving A:-Hessian operators (see [35] and the references therein). 

Just as we have for the prototypical elliptic equation, we show that 
bounded and decaying positive solutions of (HI) exhibit only two rates of 
decay: the fast decay rates and the slow decay rates. Here we say a decaying 
solution {u,v) of system (|l.ll) decays with the slow rates as |x| ^ 00 if 

u{x) ~ and v{x) ~ 

where 

„ _ /37(7-l+<?)+(7-l)o-i+o-2g __] „ _ /37(7-1 +p)+(7-1)o-2+o-iP 

Qo - - pg-{j-iy^ - Po - --■ 

On the other hand, we previously established the following equivalent char¬ 
acterization of the fast decaying solutions and whose definition is contained 
in the last statement of the theorem. 

Theorem 1 (Theorem 1 in [35]). Let 7 G (1,2], > p > 1 and ui < (12 < 0 

and let {u, v) be a positive solution of the integral system dni) where 

qo+Po< -(1-6) 

7-1 

Then the following statements are equivalent. 

(a) (n, u) € L'’°(M”) x L®°(M”), where 

n , n 

vq = — and So = —. 
qo Po 


(b) {u,v) € L''(R’") X L"(M”), where 


- 1 ) 

n — fij 


and s > max 


{ 


^(7 - 1) nj'y - 1) ^ 

n - ^7 ’ - (/37 (T2) ^ 
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(c) {u,v) is bounded and decaying, and it decays with the fast rates as 


oo, i.e., 


and 


v{x) 


u{x) ~ |x| 


-/37 


7-1 
n — /37 


T-l (In |x|) T"-! . 

p( ) -(^7+^2 ) 

7-1 


n — 13'y 

7-1 


if p(7?) 
if p(^) 
if p(=^) 


- cj 2 > n\ 

- 0-2 = n; 

- (72 < n. 


In view of this, we say a solution (tt, x) of system dni) is an integrable 
solution if {u, v) G x and say it is an optimal integrable 

solution if {u,v) € L'’(]R"') x L®(M”) for all (r, s) satisfying ()1.7p . 

There are some important observations that should be made here. The 
previous theorem asserts that solutions, under a fairly mild integrability 
assumption, are indeed bounded and fast decaying. The assumptions that 
q > p and cri < 02 are due to the in-homogeneity of the system (an issue 
which does not occur in the scalar case), but they are not so crucial. More 
precisely, the theorem still holds p > q and 02 < ui provided that these 
parameters along with u and v are interchanged in the statement of the 
theorem. Condition (ll. 6 p . which is equivalent to perhaps the more familiar 
condition 

n + ai ^ n + a 2 ^ n - 13-f 

q + 'y — l p + 'y — l~ 7 -I 

when pq > (7 — 1 )^, is certainly stronger than the condition max{( 7 o, Po} < 
(n — / 37)/(7 — 1); however, making this stronger assumption is not without 
proper motivation. For instance, in the special case of system dm) with 
a G (1,?^), it turns out that the system admits neither a bounded and 
decaying positive classical solution nor a positive integrable solution in the 
subcritical case [36] . 


n + ai 
q + l 


n + a2 

p + 1 


> n — a. 


In fact, when ai,a 2 = 0, system ()1.4p admits a positive integrable solution 
if and only if the critical case holds [ 2 ll[ 2 i|, i.e., 


1 1 n — a 

-7 ^-7 ~ - 

g -|- 1 p -|- 1 n 
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Of course, we conjecture that the same holds true for the more general Wolff 
type integral systems, but a proof of this escapes us at this time. However, 
we do have a closely related result below (see Theorem [2]). Even for system 
(HID such a result on the non-existence of positive classical solutions is 
quite non-trivial and it is often called the (generalized) HLS conjecture [3] 
(in the case of system (ll.Sp with = 1 , it is more commonly called the 
Henon-Lane-Emden conjecture). It is crucial to note that no boundedness 
or growth assumptions are being imposed here and the known partial results 
are limited to dimension n < 4 or under the aforementioned assumptions 


[MllSHlEQlEIlIM]. 

In addition, the intervals of integrability in ()1.7p are optimal. Namely, 
if a solution of (jl.lll belongs to L^(M"^) x then necessarily (see [35]) 



Therefore, q > p, cri < a 2 and (jl.Op imply that 



We also mention some previous articles concerning the analysis of the HLS 
and Wolff type integral systems, especially since they have ultimately in¬ 
spired the work in this paper. The study on the existence, non-existence, 
classification and decay properties of HLS type and related systems can be 
found in [21 jS] jS] E] (TUI [12112211221 El] • Similar studies on the Wolff type 
integral systems can be found in [H El E?! [HI [29l |23| . 

We are now ready to state the main results of this paper. 

1.1 Optimal Liouville theorem for positive solutions 

Our first main result for system (HD is the following non-existence result. 

Theorem 2. Let (3 > 0 and 7 > 1 with (5^ < n. Then the integral system 
(HD has no positive solution for any double bounded coefficients ci{x) and 
C 2 {x) if either 
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If, in particular, 7 € (1,2], then the same conclusion holds in the endpoint 
case: 

pq> {j- if and max jg'OjPo} = i ■ 

Remark. As a consequence of Theorem\^ when given some positive solution 
(u, v) to either system dni) or the quasilinear systems considered below, we 
shall therefore assume that 


Pq > {1 — if o.iT-d max|(7o,Po| 


< 


n — fl'y 

1-1 ' 


(1.9) 


The next theorem indicates that, under the assumptions in Theorem 

is indeed sharp. 


Theorem 3. Suppose that p,q,ai,a 2 , fi andj satisfy condition (II.9p . Then 
there exist double bounded coefficients ci{x) and C 2 {x) such that the integral 
system o admits a positive solution. 


1.2 Decay rates of positive solutions 

The following results concern the decay properties of bounded and decaying 
solutions starting with a result on the slow decaying solutions. 

Theorem 4. Let /3 > 0 and 7 > 1 with jl'y < n and suppose {u, v) is a 
bounded positive solution of (H]). Then the following statements hold. 

(i) If 9i < qo and 62 < Po, then there does not exist any positive constant 
c for which either 

u{x) > c(l + or v{x) > c(l + \x\)~^^ for a.e. x G M”. 

(ii) If 63 > qo, 04 > Po cmd {u,v) is not integrable, i.e., either u ^ 

or V ^ L^°(M"'), then there does not exist any positive constant C for 
which either 

u{x) < C{1 + or v{x) < C{1 + \x\)~^'^ for a.e. x G M”. 

(Hi) If{u,v) is a decaying solution but is not integrable, then it necessarily 
decays with the slow rates as \x\ — 00 , i.e., 

u{x) ~ and v(x) ~ 
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Parts (i) and (ii) of the theorem, in a sense, imply that bounded positive 
solutions that are not integrable have almost the slow rates. Then part (hi) 
demonstrates that if it is also decaying, then it actually has the slow rates. 
So this result naturally complements Theorem [H however. Theorem [5] is far 
more general than Theorem [T] in terms of the assumptions placed on system 
(HI), namely, on the parameters 7 ,p, q, ai and 02 - This leads us to ask if we 
may establish another version of Theorem [T] that relaxes the assumptions 
on the parameters. In view of this, we do have such a result, however, we 
must restrict our attention to the optimal integrable solutions. 

Theorem 5. Let /3 > 0 and 7 > 1 with fd'y < n, q> p, ai < <72 and suppose 
{u, v) is a bounded and decaying positive solution of system o satisfying 
The following statements are equivalent. 

(i) {u,v) is an optimal integrable solution. 

(ii) {u,v) decays with the fast rates as |x| —>■ 00 . 

Remark. Under the assumptions of TheoremUl the optimal integrable so¬ 
lutions are equivalent to the integrable solutions. The proof of this relies 
on key comparison estimates between the Riesz and Wolff potentials, 
the weighted HLS inequality, and a delicate bootstrap or lifting technique. 
However, we do not know if this equivalence remains true under the more 
general conditions of Theoreml^ even under the additional assumption that 
the solution is bounded and decaying. 

1.3 Quasilinear systems 

In establishing the asymptotic and Liouville type results for our family of 
integral systems, we also obtain analogous results for the general quasilinear 
system of the form 

f —divA{x,Vu) = ci{x)\xU^v‘^, 

I ( 1 . 10 ) 

y —di-vA{x,Vv) = C2{x)\x\^'^uP, 

where the map A : M” x 1 — R” satishes the following properties. The 
mapping x 7l(x,^) is measurable for all ^ G R”; the mapping ^ (-)> 7l(x,.^) 
is continuous for a.e. x G R""; for some positive constants ki < k 2 there hold 
for all ^ G R” and a.e. x G R"’, 

(a) A{xA) A > ki\Cp, 

(b) |^(x,OI < 
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(c) {A{x, — A{x, ^')) • — ?0 > 0 whenever ^ 

(d) ^(x,A^) = A|A|'^“^^(x, for all A 7 ^ 0. 

We say (n, v) is a (weak) solution of (ll.lOp if u and v belong to (M"') n 
C'(M”') and satisfy the equations in the distribution sense. In the simple 
case where A{x,^) = div^(x,Vu) is just the classical 7 -Laplace 

operator. 

To illustrate the relationship between quasilinear operators and Wolff 
potentials, we recall a consequence of the global pointwise estimates of 
[TB] [Corollary 4.13]. Namely, if {u,v) is a positive solution of (jl.lOp sat¬ 
isfying 

inf u = inf u = 0 , 

then there exist positive constants Ci and 6 * 2 , depending only on n and 7 
and the structural constants ki and ^ 2 , such that 


f C'iWi,-^(ci(y)|y|'^W'?)(x) < u{x) < C'2ld"i,^(ci(y)|y|'^W'?)(x), . . 

1 CiWi,^{c2{y)\yr^uP){x) < v{x) < C2Wi,^ic2{y)\yr^uP)ix). ^ ^ 


In view of this and by applying our results on the Wolff type integral systems, 
we shall establish the following. 


Corollary 1. Let /3 = 1 and 7 € (1,2]. For any pair of double bounded 
coefficients ci{x) and C 2 {x), system (fTTOll has no positive solution {u, v) 
satisfying 

inf u = inf u = 0 , 

R" R" 

whenever pq € ( 0 , (7 — 1 )^] or if pq > (7 — 1 )^ and 




> ’^-7 

- 7 - 1 - 


We have the following decay properties of solutions for the quasilinear 
systems. We begin with a corollary of Theorem [1] for completeness sake. 


Corollary 2 (Corollary 2 in [3S]). Let /3 = 1, 7 € (1,2], g > p > 1, 
fXi < cr 2 < 0 and let {u, v) be a positive solution of system (fTTOl) with 
Qo + Po < Then {u,v) is an integrable solution if and only if {u,v) is 

bounded, decaying and decays with the fast rates as jxj —>■ 00 , i.e., 


n — 'y 

u(x) ~ jxj WT 
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and 


if P(^) - 0-2 > n; 
if P(^) - 0-2 = n; 

if P(^) - 0-2 < n. 

Corollary 3. Lei /3 = 1, 7 € (l,n) and let {u,v) be a bounded and decaying 
positive solution of p.lUl) . If {u,v) is not integrable, then it necessarily 
decays with the slow rates as \x\ —>■ 00 , i.e., 

u{x) ~ and v{x) ~ \x\~^°. 

Corollary 4. Let /3 = 1, 7 G (l,re), q > p > I, (Xi < 02 and let {u,v) be a 
bounded and decaying positive solution of system (|1.10p satisfying Qo+po < 
Then {u,v) is an optimal integrable solution if and only if{u,v) decays 
with the fast rates as |x| —>■ 00 . 

The remaining parts of this paper are organized as follows. In ^ the 
proof of Theorem [2] is provided followed by the proof of Theorem [3l Then, 
^ and @ respectively, contains the proof of Theorem [J] and Theorem [5l 
Lastly, ^contain the proofs for the corresponding results on the quasilinear 
systems. 

Remark (Notation). Throughout this paper, we adopt the standard conven¬ 
tion that c, C, Ci,C 2 - ■ ■, are positive universal constants in the inequalities 
that may change from line to line (and sometimes within the same line it¬ 
self). 


v{x) 


n — y 

7-1 


n —7 
7 


-1 (In |x|) r-i , 


P(^yZrr')-(y-^^2) 

7-1 


2 Existence and Liouville property of solutions 


2.1 Non-existence of positive solutions 

We now prove our Liouville type theorem. 


Proof of Theorem We proceed by contradiction. That is, assume there is 
a positive solution {u, v). Let \x\ > R for some suitable R > 0 and note that 


0 < Cl < 


L 


Br{ 0) 


\yPv‘^{y)dy < C 2 < 00 . 


Then, from the first integral equation there holds 
'/sh(o) dt 


u{x) > C 


L 


|a;|+_R 


t^-0'y 


> C 


dt C 

/ t — > 1—i—■ 

l\x\+R t 
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where ao = (n —/37)/(7 —1). Inserting this into the second integral equation 
yields for |x| > R, 


v{x) >c r{[ 

J 2 \x\ ^Je 


dy 






-y'—]_ dt 

T 


> c 


paa-<T2-Pl dt 


7-1 


/2|x| 


t 


Now, if pao — ct 2 — /37 < 0, then the previous estimate implies v{x) = oo and 
we arrive at the desired contradiction. Otherwise, if pa^ — cj 2 — /I 7 > 0, we 
deduce instead 

v{x) > C\x\~^°, for \x\ > R, 

where bo = {pao — 02 — /37)/(7 — !)• Likewise, using the previous estimate, if 
qbo — cTi — /I 7 < 0 then u{x) = 00 ; otherwise, if qbo — iTi — /37 > 0, inserting 
the last estimate into the first integral equation yields 


u(x) > Clxl for |x| > R, 

where ai = (qbo — ai — (3'y)/('y — 1). Assuming that we can continue this 
procedure, we arrive at 

u{x) > C'|x|““'' and v(x) > C'|x|“^'', for |x| > R, 


where 


h 


pak -(T2- Pi 
7-1 


and Ofc 


qbk-i - ui - /?7 

7-1 


for A: = 1, 2, 3,... . 


Actually, using the definitions of and b^, we calculate that 


_ qbj_i -ai- fdj _ pq /37(7 - 1 + g) + (7 - l)cri + qa 2 

7-1 “ (7-1)2“^-' (7-1)2 

_f pq ^^2 _ 1 + g) + (^ _ l)cr^ + go- 2 , pq \ 

V (7 - 1)2 ) “^-2 (7 - 1)2 V (7 - 1)2 / 


_ 


= r-oao 


/37(7 -l + q) + ('y- l)ai + qa 2 
(7-1)2 


i=0 


0 ) 


where j = 1, 2,3,... , and 


ro 


pq 

(7-1)2- 


( 2 . 1 ) 
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We claim that this iteration process must stop after a finite number of 
steps. To see why, consider the two cases: when pq € (0, (7 — 1)^] and when 

pq> il- 1)^- 


Case 1: Suppose pq G (0, (7 — 1)^]. If pq = — 1)^, then ()2.ip implies 

that 

./37(7 - 1 + g) + (7 - l)o-i + qa2 
aj -ao j ^^2 

-00 as j — 00 . If pq G (0, (7 — 1)^), then (|2.1I) implies 


Therefore, aj, bj 
that 


aj = rguo - 
Sending j ^ 00 yields 


_ J - 1 + O') + (7 - l)f^i + g '<72 1 - rl 


(7-1)^ 


1 - ro 


< 0 and thus bj —>■ {pqo — ct 2 — / 37)/(7 — 1) < 0. 


In either case, we can find a suitably large jo such that aj^^bj^ < 0 and this 
implies u{x),v{x) = 00 , which is impossible. 

Case 2 (a): Let pq > {'j — 1 )^ and max|go,Po| > ^77^- 
Hereafter, we denote 

M = max|go,Po}- 

Let us first assume M = qo- By virtue of (12.ip . we can find a large jo such 
that 


“to = ^o°“o 


/ 37(7 - 1 + O') + (7 - l)fTi + qa 2 r^o ~ ^ 


= r; 


Jo 


n — Ij'-f 
7-1 


qo 


(7 - 1)2 

) + go < 0 . 


ro - 1 


= ri%ao - M) + M 


Thus, u{x) = 00 and we have a contradiction. Likewise, if pg > (7 — 1)^ but 
instead M = po, then we can also apply the previous iteration argument to 
deduce a contradiction. 

Case 2(b): Let pg > (7 — 1)^ and M = where 7 € (1,2]. Without 
loss of generality, we assume M = po- By Holder’s inequality, 


fj 

Jo J Bt{x) 


\yV^v'‘{y)dydt 




^ \yV^v‘i{y)dy 




t) 
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Thus, for X G i?j:j/4(0), 


l{x) >C 

Jo 




n —/S'Y+l 

> CR 


n —/S'Y+l 

> CR 


71-/37 

> CR~^i^ 


fn-P'y 

( r / \yV^v^{y)dydt)^^ 

^Jo JBt{x) ^ 

\y\'^v^{y)dydt 


|3;|+_R/4 JBt{x) 

\vV^v‘^{y)dy 


1 

7-1 


^r/a{^) 


1 

7-1 


As a result, we obtain 

71-/37 _ 

uF{x) > CR T"-! ^ 
Similarly, we can show that 

n — 13'Y ^ 

v^{x) > CR T-i ^ 


^^H/4(0) 


\yV^v'^{y)dy 


^fl/4(0) 




P 

7-1 


9 

7-1 


( 2 . 2 ) 


(2.3) 


If we multiply (12.21) by integrate over i?^/4(0)\lle(0) for suitably small 
e > 0, apply ()2.3I) then send e —^ 0, we get 


f \x\^^u^{x)dx 

J Bii/a{ 0 ) 


> 


C 


— 71 —/37„ „ I n — B'y _ ^^P 717? \ f 

7-1 P *^2 (7-1)2 ^-1 7-1 '^-Bj^/4(0) 


O f \ 

|x|°'^M^(x) dx j 


> 


O f . P9 

' \xruP(x)dx)^^, 

B„m(0) ' 


'^R/4(0) 

where the above positive constant C is independent of R since 

n — B'y n — B'y ui p np 

- -a2-n+ - - :^pq ---- 

7 — 1 (7 ~ 1) 7 “ 1 7 ~ 1 

pq — {'j — iB fn — jB^ 


-irfn-/?7 'I 

—i-7vr-«) = "■ 


7-1 

Thus, sending R ^ oo implies that \x\^'^u^{x) € L^(M”'). If we repeat the 
previous argument but instead we integrate over i?^/ 4 ( 0 )\i?jj/ 8 ( 0 ), then 


f (x) dx > C( [ \x\^'^vF{x) dx^ 

J ^R/4{0)\Bji^g{0) JBr//A(0) 


PQ 

(7-1)^ 


Br/4{0) 
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where C is independent of R. Hence, sending R ^ oo yields 



x\^'^vP{x) dx = 0. 


This implies ti = 0 and we deduce a contradiction. This completes the proof 
of the theorem. □ 


2.2 Existence of solutions 


Proof of Theorem O Indeed, we find double bounded coefficients with the 
positive radial solution pair 


u{x) 


1 

(1 + 


and v{x) 


1 

(1 + |x| 2)'^2 ’ 


where the rates 9i and 62 are specihed shortly below. In fact, for complete¬ 
ness sake, we provide a solution pair with the slow decay rates and another 
pair with the fast decay rates. 

(i) Choose the slow decay rates: 


201 = go and 202 = Po, 


so that P'y < 2p6i — a 2 < n and (5^ < 2g02 — ai < n. For \x\ < R with 
a suitable choice for i? > 0, it is obvious that u{x) and v{x), respectively, 
are proportional to Wii^ry{\y\^^v^){x) and Wii^^{\y\^^uP){x). Thus, we may 
restrict ourselves to the case where |x| is suitably large. Consider the split¬ 
ting 


Wp,,{\yrv<i){x) 



Is/Fi u 

+n—8'y ' 


H 1 +H 2 . 


Notice that for y G Bt{x), 


|a:|/2 < \y\ < 3|x|/2 whenever |x|/2 >t>0, 
and since — 0i = 0, there holds 


1 , , ,r, 9 S 2 - 

Hi > —(1 |x|^) T-i 


1 , , ,9, 9^2--^ 

> —(1 |xf) T-i 

o 


l/I f 

' i 

ill f 

' i 


M^— - 11 —h dt 
\Bt{x)\'^-^t T-i — 

1^1/2 dt \ o 2fl2g-gl-^7 1 

fT-l - > —(1 -I- |xP) 2(7-1) > ui^x). 

t o 
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Similarly, there holds 


Hi < C(l + |x| 


9g2-°'l/2 i3j 


7-1 


L 


2 fl 2 g-a '1 -/37 


tr-l — < (7(1 + |xf) 2(7-1) <Cu{x). 


Hence, we have that 


C-^Hi < uix) < CHi 


(2.4) 


for some positive constant C. On the other hand, there holds 

r 1 , 

(i+|j/| 2 ) 9«2 — dt 


ddo = 


'N/2 


fn-^'Y 


< c 


poo 

J\x\/2 




< C 


fn-yj 

oo ji+\^\ ^n+cn—26oo dr ^ ^ 

t 


l\x\/2 


fn-^'y 


r°° <ri+Pl-2qe2 dl 262 - 0 - 1-/37 

<C t r-1 — <(7(l + |xp) 2 ( 7 - 1 ) <Cu{x). 


'\x\/2 


t 


If t > |x|/2, we have that |x|/2 < \y\ < 3|(c|/2 for y G C i?j(x) and 

thus 


H2 > 


> 


>\x\/2 


I^kl/ 2 W (1+|i/| 2)'?''2 “dx dt 


fn-pi 

1 , , „ 2929-<Tl-n rOO n-/37 1 / s 

— (l + |xf) 2(7-1) / t r-1 — > u{x). 

O 4|x|/2 ^ O 


Hence, (7 ^H 2 < ■u(x) < CH 2 for some positive constant (7, and by com¬ 
bining this with (j2.4|] . we obtain 

u(x) = ci(x)IHa,7(|dr^^^'^)(a;) 

for some double bounded function ci{x). Likewise, similar calculations on 
the second integral equation will lead to 

v{x) = C 2 {x)W^^^{\ypvF){x) 

for some double bounded function C 2 {x). 
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(ii) Choose the fast decay rates: if the stronger condition p > 

and q > hold, then we can take the fast rate 26i = 202 = 

in which 20ip — (T 2 > n and 202 ? — (Ji> n. Then 

Wf}^^{\y['^v‘i){x) = - = ci{x)u{x), 

(1 + |x|2) 2(7-1) 

Wp,^{\yV^vP){x) = - = C2{x)v{x), 

(1 + |xp) 2(7-1) 

for some double bounded functions ci(x) and C 2 (x). This follows from direct 
calculations as before and so we omit the details. This completes the proof 
of the theorem. □ 


3 Slow decay of positive solutions 


In this section, (u, v) is always taken to be a bounded positive solution of 
system (HI, and the goal is to prove Theorem SI 

Proposition 1. Let 0i < and 02 < Po- Then there does not exist any 
positive constant c for which 

either u{x) > c(l + or v{x) > c(l + for a.e. x € M"'. 


Proof. The proof incorporates the iteration scheme found in the proof of 
Theorem [21 On the contrary, assume there exists a positive constant c such 


that 


u{x) > c(l + |x|) 

Indeed, for large |x| and with (0), there holds 

\yr{l + \y\)-<>^Pdy^^dt 


v{x) >c r ( [ 

J2\x\ 


in-p'y 


> c(I + |x|) 


-ai 


where 0o = and oi = Ml_£ 2 _ 22 .^ Inserting this estimate into the first 
integral equation yields 


u{x) > c 


r f 

|yrHi + |y|) 

/2\x\ ^ 

fn-p-t ) 


)-'->*>e(l + |x|r*>, 


where bi = As before, provided that the rates remain positive, 

we can repeat this procedure inductively to arrive at 


v{x) > c(l + |x|) and u{x) > c(l + |x|) 
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where 


phj - /37 - 0-2 qaj - /37 - fii 

Oj+i = — - --- and bj = 


for j = 1,2,3,... . 


7 — 1 7 ~ 1 

If we set ro = and 70 = /37(7 “ 1 + ?) + (7 “ l)^i + ^^29; we calculate 

that 




_ Qaj - /37 - fJi _ 1 r pbj_i - fi'y - a 2 


7-1 

pqbj-i - 70 ^ 

(7 - 1 )^ (7 - 1 ) 


i{ 


7-1 

qaj-i - P'y -ai 

7 - 1 


P'y - o-ij 


70 


= robj -2 - 


Vo , N i, Vo 

(1 + ro) = ... = r^5o - 


(7 - 1) 


(7-1)2" 


i=0 


= - VO - 7 -TTo = ri{9i - qo) + go- 

pq-i'j- ly 

Since 6*1 < gO) this implies that we can choose a suitably large jo such that 
bjQ < 0. Therefore, 


roo dt 

v(x) >C t r-l — = 00 , 

^ ^ “ J 2 \x\ t 

but this is impossible. Likewise, if there exists a positive constant c such 
that v{x) > c(l + |x|)“^ 2 , then we can apply the same iteration scheme to 
deduce a contradiction. This completes the proof of the proposition. □ 


Proposition 2. Let 63 > go and 6*4 > po and {u, v) is not integrable. Then 
there does not exist any positive constant C for which 

either u{x) < (7(1 + or v{x) < (7(1 + \x\)~^^ for a.e. x G M”. 

Proof. Assume there exists a (7 > 0 such that u{x) < (7(1 + |x|)“^®. Then 
n — (n/go )03 < 0 and for a suitable choice of i? > 0 , we obtain 


u{xY° dx = / u{x)‘>odx= / u{x)‘>odx + 
Jm.^ Jbr{0 ) 

< (7i + C 2 [ — < 00 . 

Jr t 




u{x) 'JO dx 


Similarly, if there exists a (7 > 0 such that v(x) < (7(1 + 1x1)“^"^, then 
n — (u/po )^4 < 0 and we can show v € L^°(M"'). In any case, we arrive at a 
contradiction with (u, v) being not integrable. □ 
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Proposition 3. If {u, v) is not integrahle but is a decaying solution, then 
{u,v) necessarily decays with the slow rates as \x\ —> oo. 

Proof. This follows immediately from Propositions [T] and [ 2 j □ 

Proof of Theorem The theorem follows from Propositions HHSl . □ 


4 Fast decay of positive solutions 

This section contains the proof of Theorem [5] and throughout the section, 
we assume the same conditions as those stated in the theorem. Moreover, 
(u, v) is assumed to be an optimal integrable solution of unless stated 
otherwise. 


Proposition 4. For suitably large |x|, there exists a constant c > 0 such 
that 


u{x),v{x) > c\x\ 


71 — 

7-1 . 


Proof. For suitably large \x\, we have that 

f 

u{x) >c 

Jl 


/ bi ( o ) \yf^^'^iy)dy ^^dt ^ ^ 

i+ixi ^ ^ t - 


°° P~t-n (If I I — 

t T"-! - > C\x\ T"-! . 


The corresponding estimate for v{x) follows similarly. 


□ 


Proposition 5. There holds u{x) ~ |x| 


7-1 


Proof. By the hypotheses, there is some rate 6 i for which u{x) ~ |x| 
First, we claim that 9i > ; otherwise, we can find e > 0 so that 9i < 


If we set r = , then for a sufficiently large i? > 0, 


■>) > c 


/ 

JK’ 


\Br(0) 


Jr 


t 


= oo. 


This is impossible since dni) ensures that u G U 
claim. Then Proposition [J] further implies that 


), and this proves the 


01 = 


n — 13'y 

7-1 ■ 


□ 

^ _ n — 0'y 

Proposition 6 . If — a 2 > n, then v(x) ~ \x\ t"-! . 
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Proof. As before, there is some positive rate 62 such that v{x) ~ \x\ By 
the integrability of u, we claim that 


■ n-/37 - (/37 + ^2) -| 

7 — 1 ’ 7 “ 1 


I- 


62 > min|- 
However, notice that 

- /37 - (/37 + 0-2) 'i n - /37 


min 


jn| 


n 


7-1 


7-1 


7 - 1 ’ 


since — (72 > n. Assume the contrary. Then we can find e > 0 such 


that 


^ / ra - /?7 n (7 - 1 + e) 

62 < -:- and set s = 


7 — 1 + e n — j3'y 

We can choose i? > 0 suitably large so that 


IS 

\l^ 


^)>C 


/ 

Jw 


dx> c 


\Br{0) 


t n-Pl - = CXD, 


'R 


but this is impossible since (|1.7p ensures that v G This proves the 

claim and by combining this with Proposition 01 we obtain that 


62 = 


n — 13'y 

7 - 1 ■ 


□ 


Proposition 7. If — 02 = n, then 

Proof. We shall make use of the following identity which follows from ele¬ 
mentary arguments from calculus (see m)- For A > 0, 


v[x 


_ n — jS'y 1 

\x\ (In |x|) ^-1 , 


n — (3'y 


lim 


x\ T"-! 1"°° ^ Inf dt 7 — 1 

t n — (3^ 


(lnA|x|)T'-i Jm^\ 
For A G (1/2,1), we write 


(4.1) 


v{x) < C 



\yV^uP{y)dy 



n — iB'y 

For large |x|, since u{x) ~ |x| and 

n — (3^ a 2 + (3 j n — (3^ 

^(7 — !)^ 7“1 7~1 ' 


Cih+h). 
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we have that 


_ n-l3~f . <T 2 fM^\ fjf „ n-0'i , <T2 +^j n-fl^ 


h < C\x\ 
Thus, 


n — 13'y _ 1 

lim |x| T-i (ln|x|) t-i/i = 0. 

|a;|—>oo 


(4.2) 


Likewise, choose a sufficiently large R > 0. Then for large |x| we have that 


h<c 


< c 


< c 


< c 


fn-P'r 


°° // b «( 0 ) + \y\^^uP{y)dy 


L 
L 

f°° f Cl + C 2 r 

J\ 

j 


t 


l\\x\ 

roo 


t^-d'Y 

*+hl 1 


■y-^ dt 


fn-13'y 

Cl+ C2r-^ hr ^ ^ 

- \ (_/ 


fn-fS'i 


poo 

7 

Aja^l 


In t \ dt 


yn-M J I ■ 


Combining this estimate with identity (I4.ip and sending A —>■ 1 yields 


1^1“^°° (In |x|) T"-! 
Hence, (j4.2p and (14.3p imply 

1-/37 


n — ^'y 

lim h < C. (4.3) 


X, 7~1 

lim - —v{x) < C. 


hl^CXD |x|) 7-1 


(4.4) 


Notice that for any A > 1, i?^_|3.|(0) C Bt{x) if f > A|x|. For a proper choice 
of i? > 0, Proposition m implies that 


fiy-Pi 


, r 7^*-n(o)\bh(o) dt 

;(x) > c / ' ' 

hA|; 

>_.f 

J\\: 


l\\x\ 


> c 


poo 
J\\x\ 


In t \ dt 


tn-0-f J f 
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Thus, applying identity (|4.ip to this then sending A —1 yields 


n —/S'Y 

X 

lim - —v(x) > c > 0. (4.5) 

|x|^00 |x|) T-l 

Hence, (14.41) and (14.5p imply the desired result. □ 


Proposition 8. If p{ "' — a 2 < n, then 


)-(/3T'+g'2) 


v[x 


\x\ 


7-1 


Proof. Fix some i? > 0 and for |x| > 2R, write 

v{x) < c(^ J 


+ 


/Of 




f ^-01 


where Hi = [|x| — R, |x| + i?]. Then we claim that 

p(4^=^)-(/3t+ct2) 

(a) lim| 2 ,|^oo |2;| ~ Ji = 0, 

)-{Pr+<y2) 

(b) lim|j,|^oo |a;| J 2 = C, 

and the desired result follows once we prove this. 

(a) Indeed, Bt{x) C i?2i+_R(0) if t € Hi and so Proposition [5] yields 


Ji<C 


/ 

Jn 


Jo 


7-1 


^)dr 
r 


in-^'Y 


^ < c\x\ 
t ~ ' 


p(4^7^)-(^T+CT2) 


7-1 


-1 


and part (a) follows accordingly. 

(b) To show this part, we first verify that 



isde) 


0-2-p( 


n —/37 

7-1 


dy\T^ dt 


fn-l3y 


< 00, 


(4.6) 


where e is any unit vector. To show this, let c G (0,1) and consider the 
splitting 


J 3 + Ja 




fstie) \y\ 


0-2-p( 


n—(3‘~^ 

7-1 


dy\^ dt 


fn-fS'Y 


t 
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Indeed, 


_ / n — I3'y \ 

, ^ r ffB,ie)\y\ ~ dt ^ ^ r/\B,{e)\^^ 

* - Jo y -- ) T-'^Jo 

rc 


'y — dt 
t 


Jj_ dt 

< / tT^-1 — < oo, 


Jo t 

since y G Bt{e) ensures that l — c<\y\<l + c. We can find a suitably large 
i? > 0 so that 


Ja<C 




in-13'1 


ri?t n+cr2-p(^l^)^ _i_ 

Jo r \ r-i 


fn-p'y 
oo P{ ) - (g7+g'2) 


7-1 


rf> \ 7—1 dt 

~) T 

dt 

— < oo. 

t 


This completes the proof of the claim. Now we turn our attention to esti¬ 
mating the term J 2 - Indeed, by setting 

112 = [1 - R/\x\^ 1 + R/\x\]-, 

using the change of variables 


z = — s = — 
1x1’ 1x1’ 


and applying (14.61) . we get that 

Jo<cf 

“ V J t 


<c j 

Jq,9. 


Js4x/|a:|) 1^1 as 


gU-P'l 


■ X 


7-1 


< Cl|x 




7-1 


Likewise, using Proposition [4] combined with similar arguments as above, 
we can also show that 


J 2 > c|x 


P( 1y_^L )-(/3T+g2) 


7-1 


This 


proves part (b) and thus completes the proof of the proposition. 


□ 
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Proposition 9. Suppose that {u, v) is a bounded and decaying solution of 
system (HH). If {u,v) decays with the fast rates as \x\ oo, then it is 
an optimal integrable solution, i.e., {u,v) G x L^(R"’) for all {r,s) 

satisfying condition m- 


Proof. Suppose that {u,v) decays with the fast rates as |x| —>■ oo and (r, s) 
satisfies (HZl). 

(i) If u{x) decays with the rate , then we can find a suitably large 
i? > 0 so that 

f u{xY dx < f u{xY dx + f u{xY dx 

Jk" Jbr(O) Jr^\Br(0) 


>Br{0) 

< Cl+ C 2 


\Br{0) 


^\Br{0) 


'L 


dt 


x| dx < C ! t T--! — <00. 


t 


Likewise, if v{x) decays with the rate we can use similar argu¬ 

ments to show V € 

_ n — jS'y 1 

(ii) Assume v{x) ~ \x\ r-i (In IxDt'-i . For small e > 0, we can find R > 0 
such that 


(InlxDT-i < C\xY for |x| > R. 


Thus 


[ v{xY dx<Ci + C2 [ — < 00, 

JR" Jr I 

since n — s '^Y-i + e < 0 provided that e is chosen to be small enough, 
(iii) Assume 


n-P-i , cr2+P-y 

v{x) ~ |x| T"-! . 


Condition (II.7p ensures that 

s f n — Y'j 


n — 


Y) - (/37 + f^2)) < 0, 


and thus 


7-IV 7- 

/ u(x)'* dx<Ci + C2 [ 
JR" Jr 


dt 


< 00. 


In any case, {u,v) G L^(]R"') x and this completes the proof. 

□ 

Proof of Theorem O This follows directly from Propositions [5] to [9l □ 
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5 Systems of quasilinear differential equations 

Proof of Corollary m On the contrary, assume that {u,v) is a positive solu¬ 
tion of system (ll.lOp satisfying 

inf u = inf v = 0. 

Since the coefficients are double bounded, the global estimate (ll.lip ensures 
there is a positive constant C such that 

C-^Wi,^{\yrv^){x) < u{x) < CWi^^{\yrv^){x), 
C-^Wi^^{\ypvF){x) < v{x) < CWi^^{\ypuP){x). 

From this, we clearly get two double bounded functions ci(x) and C 2 (x) such 
that 


u{x) = ci(x)PFi,^(|y|'^iu'')(x), 
v{x) = C 2 {x)Wi^^{\y\^^vF){x), 

but this is impossible in view of Theorem [2j □ 

Proofs of Corollaries 0 an d\^ Indeed, by virtue of estimate (jl.lip . there 
exist double bounded coefficients ci(x) and C 2 (x) such that (n, v) is a positive 
solution of the integral system 

u{x) = ci(x)ITi,^(|y|'"ix'?)(x), 
v{x) = C 2 {x)Wi^^{\yy^uF){x). 

Then the results follow directly from Theorem 0] and Theorem [5l □ 

Acknowledgment: The author wishes to thank the anonymous referees 
of this paper. Their comments and suggestions have certainly improved 
the presentation and overall quality of this manuscript. Part of this work 
was completed while the author was visiting the School of Mathematical 
Sciences at the University of Science and Technology of China in Hefei, PRC 
and the Department of Applied Mathematics at Northwestern Polytechnical 
University in Xi’an, PRC. The author wishes to thank both institutions for 
their support and hospitality. 


24 








References 

[1] L. Caffarelli, B. Gidas, and J. Spruck. Asymptotic symmetry and local 
behavior of semilinear elliptic equations with critical Sobolev growth. 
Comm. Pure Appl. Math., 42 (3) (1989) 271-297. 

[2] L. Caffarelli, R. Kohn, and L. Nirenberg. First order interpolation 
inequalities with weights. Compos. Math., 53 (3) (1984) 259-275. 

[3] G. Caristi, L. D’Ambrosio, and E. Mitidieri. Representation formu¬ 
lae for solutions to some classes of higher order systems and related 
Liouville theorems. Milan J. Math., 76 (1) (2008) 27-67. 

[4] H. Chen and Z. Lii. The properties of positive solutions to an integral 
system involving Wolff potential. Discrete Contin. Dyn. Syst., 34 (5) 
(2014) 1879-1904. 

[5] Z. Cheng and C. Li. Shooting method with sign-changing nonlinearity. 
Nonlinear Anal., 114 (2015) 2-12. 

[6] W. Chen and C. Li. Radial symmetry of solutions for some integral 
systems of Wolff type. Discrete Contin. Dyn. Syst., 30 (4) (2011) 1083- 
1093. 

[7] W. Chen and C. Li. Super polyharmonic property of solutions for PDE 
systems and its applications. Commun. Pure Appl. Anal., 12 (6) (2013) 
2497-2514. 

[8] W. Chen, C. Li, and B. Ou. Classification of solutions for a system 
of integral equations. Comm. Partial Differential Equations, 30 (2005) 
59-65. 

[9] W. Chen, C. Li, and B. Ou. Classihcation of solutions for an integral 
equation. Comm. Pure Appl. Math., 59 (3) (2006) 330-343. 

[10] L. D’Ambrosio and E. Mitidieri. Hardy-Littlewood-Sobolev systems 
and related Liouville theorems. Discrete Contin. Dyn. Syst. Ser. S, 7 
(4) (2014) 653-671. 

[11] J. Davila, M. del Pino, M. Musso and J. Wei. Fast and slow decay 
solutions for supercritical elliptic problems in exterior domains. Calc. 
Var. Partial Differential Equations, 32 (4) (2008) 453-480. 


25 



[12] J. Dou and H. Zhou. Liouville theorems for fractional H&on equation 
and system in M”. Commun. Pure Appl. Anal., 14 (5) (2015) 1915-1927. 

[13] B. Gidas and J. Spruck. A priori bounds for positive solutions of non¬ 
linear elliptic equations. Comm. Partial Differential Equations., 6 (8) 
(1981) 883-901. 

[14] B. Gidas and J. Spruck. Global and local behavior of positive solutions 
of nonlinear elliptic equations. Comm. Pure Appl. Math., 34 (4) (1981) 
525-598. 

[15] M. Hmon. Numerical experiments on the stability of spherical stellar 
systems. Astronom. Astrophys., 24 (1973) 229-238. 

[16] T. Kilpelainen and J. Maly. The Wiener test and potential estimates 
for quasilinear elliptic equations. Acta Math., 172 (1) (1994) 137-161. 

[17] Y. Lei. Decay rates of solutions of an integral system of Wolff type. 
Potential Anal., 35 (4) (2011) 387-402. 

[18] Y. Lei. Asymptotic properties of positive solutions of the Hardy-Sobolev 
type equations. J. Differential Equations, 254 (4) (2013) 1774-1799. 

[19] Y. Lei and C. Li. Integrability and asymptotics of positive solutions of a 
7 -Laplace system. J. Differential Equations, 252 (3) (2012) 2739-2758. 

[20] Y. Lei and C. Li. Decay properties of the Hardy-Littlewood-Sobolev 
systems of the Lane-Emden type (2013). http://arXiv:1302.5567. 

[21] Y. Lei and C. Li. Sharp criteria of Liouville type for some nonlinear 
systems. Diserete Contin. Dyn. Syst, 36 (6) (2016) 3277-3315. 

[22] Y. Lei, C. Li, and C. Ma. Asymptotic radial symmetry and growth 
estimates of positive solutions to weighted Hardy-Littlewood-Sobolev 
system of integral equations. Cale. Var. Partial Differential Equations, 
45 (1-2) (2012) 43-61. 

[23] Y. Li. Asymptotic behavior of positive solutions of equation Aw -|- 
K{x)uP = 0 in J. Differential Equations, 95 (2) (1992) 304-330. 

[24] E. Lieb. Sharp constants in the Hardy-Littlewood-Sobolev and related 
inequalities. Ann. of Math., 118 (2) (1983) 349-374. 


26 



[25] G. Lu and J. Zhu. Symmetry and regularity of extremals of an integral 
equation related to the Hardy-Sobolev inequality. Calc. Var. Partial 
Differential Equations, 42 (3-4) (2011) 563-577. 

[26] E. Mitidieri. Nonexistence of positive solutions of semilinear elliptic 
systems in . Differ. Integral Equations, 9 (3) (1996) 465-480. 

[27] E. Mitidieri and S. I. Pokhozhaev. A priori estimates and the absence 
of solutions of nonlinear partial differential equations and inequalities. 
Tr. Mat. Inst. Steklova, 234 (2001) 1-384. 

[28] Q. H. Phan. Liouville-type theorems and bounds of solutions for Hardy- 
Henon systems. Adv. Differential Equations, 17 (7-8) (2012) 605-634. 

[29] N. C. Phuc and L E. Verbitsky. Quasilinear and Hessian equations of 
Lane-Emden type. Ann. of Math., 168 (3) (2008) 859-914. 

[30] P. Polacik, P. Quittner, and P. Souplet. Singularity and decay esti¬ 
mates in superlinear problems via Liouville-type theorems, I: Elliptic 
equations and systems. Duke Math. J., 139 (3) (2007) 555-579. 

[31] J. Serrin and H. Zou. Non-existence of positive solutions of Lane-Emden 
systems. Differ. Integral Equations, 9 (4) (1996) 635-653. 

[32] P. Souplet. The proof of the Lane-Emden conjecture in four space 
dimensions. Adv. Math., 221 (5) (2009) 1409-1427. 

[33] S. Sun and Y. Lei. East decay estimates for integrable solutions of the 
Lane-Emden type integral systems involving the Wolff potentials. J. 
Eunct. Anal, 263 (12) (2012) 3857-3882. 

[34] J. Villavert. Shooting with degree theory: Analysis of some weighted 
poly-harmonic systems. J. Differential Equations, 257 (4) (2014) 1148- 
1167. 

[35] J. Villavert. A characterization of fast decaying solutions for quasilinear 
and Wolff type systems with singular coefficients. J. Math. Anal. AppL, 
424 (2) (2015) 1348-1373. 

[36] J. Villavert. Qualitative properties of solutions for an integral system 
related to the Hardy-Sobolev inequality. J. Differential Equations, 258 
(5) (2015) 1685-1714. 


27 



[37] J. Villavert. Sharp existence criteria for positive solutions of Hardy- 
Sobolev type systems. Commun. Pure Appl. Anal, 14 (2) (2015) 493- 
515. 

[38] J. Yang. Fractional Sobolev-Hardy inequality in Nonlinear Anal, 
119 (2015) 179-185. 


28 



